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Abstract 

Background: In this researcli, tlie rennoval of natural organic nnatter fronn aqueous solutions using advanced 
oxidation processes (UV/H2O2) was evaluated. Therefore, the response surface methodology and Box-Behnken 
design matrix were employed to design the experiments and to determine the optimal conditions. The effects of 
various parameters such as initial concentration of H2O2 (100-180 mg/L), pH (3-1 1), time (10-30 min) and initial 
total organic carbon (TOC) concentration (4-10 mg/L) were studied. 

Results: Analysis of variance (ANOVA), revealed a good agreement between experimental data and proposed 
quadratic polynomial model (R^ = 0.98). Experimental results showed that with increasing H2O2 concentration, time 
and decreasing in initial TOC concentration, TOC removal efficiency was increased. Neutral and nearly acidic pH 
values also improved the TOC removal. Accordingly, the TOC removal efficiency of 78.02% in terms of the 
independent variables including H2O2 concentration (100 mg/L), pH (6.12), time (22.42 min) and initial TOC 
concentration (4 mg/L) were optimized. Further confirmation tests under optimal conditions showed a 76.50% of 
TOC removal and confirmed that the model is accordance with the experiments. In addition TOC removal for 
natural water based on response surface methodology optimum condition was 62.15%. 

Conclusions: This study showed that response surface methodology based on Box-Behnken method is a useful 
tool for optimizing the operating parameters for TOC removal using UV/H2O2 process. 
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Background problems such as biological growth, taste, odor, corrosion 

Natural organic matters (NOMs) are complex compounds and membrane fouling [4-6]. NOMs can be removed via 

with different molecular composition and size that can be several methods such as coagulation, conventional filtra- 

found in surface and ground water resources. Their quan- tion through different media, enhanced coagulation, mem- 

titles, specifications and degradation depend on the sea- branes, ion exchange,adsorption and advanced oxidation 

son, weather, microorganisms activities, human activities process [5,7,8]. While the removal of NOM in order to re- 

and the eutrophication status of water resources [1-3]. duce the formation of disinfection byproducts is becoming 

The presence of NOM in water treatment processes and increasingly important, some conventional treatment 

distribution systems is undesirable and cause several methods cannot remove NOMs effectively [7,9]. Coagula- 
tion process (including the use of iron salts and alum) is 
effective to remove high molecular weight fraction of 
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widespread in developing countries [5,7]. Advanced oxida- 
tion processes (AOPs) are a variety of methods to elimin- 
ate or oxidize NOM from raw waters. These process 
include O3/H2O2, O3/UV, UV/H2O2, Ti02/UV, H2O2/ 
catalyst, Fenton and photo-Fenton processes as well as 
ultrasound [5,10,11]. In general, advanced oxidation pro- 
cesses, refers to the processes hydroxyl radicals (OH-) are 
generated as an intermediate via different methods [12]. 
Due to highly oxidizing potential of hydroxyl radicals 
(2.8 v), advanced oxidation processes are often based on 
production of this active radical. Application of UV radi- 
ation and a suitable oxidizing agent such as H2O2 is one of 
the most effective methods for advanced oxidation. This 
method is highly effective for removing NOM and refrac- 
tory organic compounds from water. In this process, 
H2O2 molecule is divided into two hydroxyl radicals be- 
cause of absorbing photons. These radicals can attack or- 
ganic molecules under proper operating conditions to 
produce end products such as water, CO2 and inorganic 
acids [5,11,13]. During treatment with UV/H2O2, NOM 
compounds are mostly oxidized and aromatic contents 
are reduced. Also during the process, high molecular 
NOM are transformed into low biodegradable compounds 
such as aldehydes, carboxylic acids [1,14]. Generated hy- 
droxyl radicals, reduces both total organic carbon (TOC) 
and disinfection by product formation potential (DBPFP) 
in raw water [1]. It has been reported that at proper con- 
centration of hydrogen peroxide and appropriate UV dose, 
this process can completely mineralize NOM into inor- 
ganic compounds [14]. 

In most related studies, experiments have been con- 
ducted with changing the independent variables and keep- 
ing the others constant. In such studies, in addition to 
large number of runs, spending time and cost, the interac- 
tions between the operating variables are not revealed. 
These restrictions can be eliminate using the response sur- 
face methodology (RSM) [15,16]. RSM is a collection of 
mathematical and statistical techniques applied for devel- 
oping, improving and optimizing a response affected by a 
number of process variables. This method has advantages 
such as fast and reliable data achievement, understandable 
interaction effects of different parameters and conse- 
quently a significant reduction in experimental runs, 
time and cost [15-17]. RSM have been successfully ap- 
plied in different process optimization such as chemistry, 
biochemistry, physics, environmental science, membrane 
process and AOPs [15,16]. 

The aim of present work was to optimize UV/H2O2 
operating parameters for natural organic matter deg- 
radation from aqueous solutions by using response 
surface methodology. The effect of main variables in- 
cluding pH, time, H2O2 and initial TOC concentration 
were evaluated using experimental data and a predicted 
model was developed. Finally, based on RSM optimum 



condition, TOC removal was studied for natural water 
samples. 

Materials and methods 

Chemicals and UV/H2O2 system 

In this study, humic acid ( Sigma- Aldrich, Germany) and 
Hydrogen peroxide 30% (Merck, Germany) were used. 
Sodium hydroxide (Merck, Germany) and chloridric acid 
(Merck, Germany) were used for pH modification. UV/ 
H2O2 system consisted of a batch reactor made of smooth 
stainless steel with an effective volume of 2.5 L, 92 cm 
length and 7.6 cm diameter. A 55-watt low-pressure mer- 
cury lamp, ultraviolet radiation source (UV-C) with an in- 
tensity of 50 mW/cm^ (Arda company, France) was used. 

Preparation of synthetic and raw water 

In this research, TOC removal from synthetic and raw 
water samples were studied. Synthetic samples were pre- 
pared by dissolving a known amount of humic acid in 
deionized water [18]. Raw water samples were taken in 
July (2012) from Sanandaj water treatment plant (Kurdistan 
Province, West of Iran). Samples were immediately stored 
at 4°C prior to use. The average alkalinity of raw water 
samples was 200 mg/L as CaCOs, the pH was 8.15, the 
turbidity was 3.5 NTU, the TOC was 4.2 mg/L, the dis- 
solved organic carbon was 3.3 mg/L and the absorption 
coefficient of UV254 was 0.098 cm'^. 

Batch experimental program & analytical methods 

In this study, TOC was measured as the indicator for 
humic substances [1]. In each experiment, 2.5 L of feed 
solution with known initial concentration of TOC, H2O2 
and pH value was poured in the batch reactor for UV ir- 
radiation. Sampling was conducted at proper intervals 
for TOC measurement. TOC was analyzed using a TOC 
analyzer (TOC-VCPH, Shimadzu, Japan). The percent- 
age removal (Y) of TOC was calculated by Equation 1. 

Y (%) = 100 (Co-Cf)/Co (1) 

where, Co is the initial concentration of TOC (mg/L), Cf 
is the final concentration of TOC (mg/L) after a known 
time. 

Experimental design and data analysis 

For determination of independent variables affecting 
the TOC removal by UV/H2O2 process and their inter- 
action effects, RSM method based on Box-Behnken design 
was applied [17,19,20]. For this design, the main variables 
are H2O2 initial concentrations (100-180 mg/L), pH 
(3-11), time (10-30 min) and initial TOC concentra- 
tions (4-10 mg/L). Independent variables in the range 
of -1 and +1 were coded according to classical methods 
similar to previous studies (Table 1). In RSM for each 



Rezaee et al. Journal of Environmental Health Science & Engineering 2014, 12:67 
http://www.ijehse.conn/content/1 2/1/67 



Page 3 of 8 



Table 1 Independent variables and experimental levels 



Variable 


Factor 


Unit 




Range and level 








Low(-I) Middle (0) High (+1) 


H2O2 concentration 


Xi 


mg/L 


100 


140 180 


pH 


X2 




3 


7 11 


Time 


X3 


min 


10 


20 30 


Initial TOC 


X4 


mg/L 


4 


7 10 



response a model is defined that can predict the individ- 
ual and interaction effect of different parameters. Gen- 
eral form of a quadratic model for four variables is as 
Equation 2 [19,21]. 

= P° + |3iXi + 132X2 + 133X3 + 134X4 + |3iiXi2 + 1322X2^ + 1333X32 + 1344X4^ 

+ Pi2^1^2 + P13X1X3 + |3]^4XiX4 + 1323X2X3 + 1324X2X4 + 1334X3X4 

(2) 

Where Y = predicted response 
|3o = constant coefficient 
Pi> p2> Ps and P4 = linear effect coefficients 
Pii> p22> p33> and P44 = quadratic effect coefficients 
Pi2>Pi3> Pi4> P23> P24 and, P34 = interaction effect coefficients 
Xi, X2, X3, and X4 = independent variables 
Design-expert (Stat-Ease, trial version) was the software 
used for designing the experiment, statistical analysis 
(e.g. analysis of variance) and response surface studies. All 
presented graphs were generated using the software. 

Results and discussion 

Development and evolution of prediction model 

The experimental results, the model predictions and the 
response (Y) using Box-Behnken design matrix for UV/ 
H2O2 process is presented in Table 2. The results of 
analysis of variance (ANOVA) for the independent vari- 
ables are presented in Table 3. The quadratic polyno- 
mial model (Equation 3) for the dependent variable (Y) 
was fitted. Based on the RSM and quadratic polynomial 
equation, there can be an empirical relationship be- 
tween the response and the independent variables. The 
final obtained equation based on the coded factors is as 
follows. 

Y = 71.50 + 3.56x1-15.08x2 + 3.58x3-11.19x4 
-0.75xiX2 + 2.25x1X3 + 0.075x1X4 + 4.00x2X4 
-1.50x3X4 + 0.52x1^-16.51 X2^-2.76 xs^-SA? x^^ 

(3) 

Routinely to assess the adequacy of a model, the coef- 
ficient of determination (R^) and the lack of fit test is 
used [22,23]. Coefficient of determination (R^) refers to 
the changes described by the model to the overall 
changes. Therefore, whatever R^ is closer to 1, the power 



Table 2 The Box-Behnken design matrix for variables 
along with Experimental and predicted values response 

Run Xi X2 X3 X4 TOC removal efficiency (Y) 



Experimental (%) Predicted (%) 



1 


180 


7 


20 


4 


84 


82 


2 


180 


7 


30 


7 


80 


78.65 


3 


140 


3 


20 


4 


78 


79.79 


4 


140 


7 


10 


4 


70 


69.73 


5 


140 


7 


30 


10 


53 


54.15 


6 


100 


7 


30 


7 


69 


67.04 


7 


140 


3 


20 


10 


52 


49.40 


8 


140 


7 


30 


4 


82 


79.54 


9 


140 


7 


10 


10 


47 


49.99 


10 


140 


11 


10 


7 


35 


33.56 


11 


180 


7 


20 


10 


61 


58.99 


12 


140 


7 


20 


7 


70 


71.50 


13 


140 


7 


20 


7 


71.5 


71.50 


14 


100 


7 


20 


10 


48.5 


51.73 


15 


100 


7 


10 


7 


64 


64.37 


16 


100 


7 


20 


4 


71.8 


74.26 


17 


100 


11 


20 


7 


40 


37.62 


18 


140 


11 


30 


7 


38 


40.73 


19 


100 


3 


20 


7 


68 


66.29 


20 


140 


3 


10 


7 


66 


63.73 


21 


180 


3 


20 


7 


72 


74 


22 


140 


3 


30 


7 


69 


70.89 


23 


180 


11 


20 


7 


41 


43.24 


24 


180 


7 


10 


7 


66 


66.99 


25 


140 


11 


20 


4 


40 


41.62 


26 


140 


11 


20 


10 


30 


27.24 


27 


140 


7 


20 


7 


73 


71.50 



of fitted model is greater to describe the response 
changes as a function of the independent variables [24] . 

Based on the ANOVA results, R^ for the removal of 
TOC is 0.98 that shows the adequacy of built model. 
Lack of fit test is a sign of lack of experimental data for 
a model that the model can not calculate the random er- 
rors of experimental data. If the test is statistically sig- 
nificant (p value < 0.005), it indicates that the model is 
not fit for the prediction of response [22,23]. With re- 
gard to Table 3, insignificant value of lack of fit test 
(p value > 0.05) relative to the pure error indicates that 
there is good correlation between the variables and 
process response. Besides determination of coefficient 
(R^), adjusted determination of coefficient (R^ adjusted) is 
a useful statistical tool to evaluate the model adequacy 
[15,22,23]. Radjusted is high (0.96) and close to R^ value, 
indicates the adequacy of the developed model to 
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Table 3 Analysis of variance (ANOVA) results of quadratic model for TOC removal efficiency (Y) 



Source 


Sum of squares 


Df 


Square mean 


F value 


p-value prob > F 


Remark 


Model 


6336.57 


14 


452.61 


46.56 


< 0.0001 


Significant 


Xi 


151.94 




151.94 


15.63 


0.0019 




X2 


2730.08 




2730.08 


280.84 


< 0.0001 




X3 


154.08 




154.08 


15.85 


0.0018 




X4 


1503.04 




1503.04 


154.62 


< 0.0001 




Xi X2 


2.25 




2.25 


0.23 


0.6391 




X] X3 


20.25 




20.25 


2.08 


0.1745 






0.023 




0.023 


2.315E-003 


0.9624 




X2 X3 


0.000 




0.000 


0.000 


1 .0000 




X2 X4 


64.00 




64.00 


6.58 


0.0247 






9.00 




9.00 


0.93 


0.3549 




X? 


1.47 




1.47 


0.15 


0.7042 




A2 














xi 


40.70 




40.70 


4.19 


0.0633 




y-i 


159.87 




159.87 


16.45 


0.0016 




Residual 


116.65 


12 


9.72 








Lack of fit 


112.15 


10 


11.22 


4.98 


0.1786 


Not significant 


Pure error 


4.50 


2 


2.25 








Cor total 


6453.22 


26 










Note: R^ = 0.9819; Adj 


= 0.9608; C.V% = 


5.13 and PRESS 


= 656.13. 









predict the process response (Table 3). According to 
Table 3, for response (Y), the quadratic model was sta- 
tistically significant (P <0.0001). Other significant terms 
of the model are X„ X2, X3, X4, X2X4, xi xi Normal 
probability plot is a graphical method for determining 
residuals normality [21,25]. A normal probability plot of 
the residuals versus the response (Y) is presented in 
Figure 1. Graphical data on the plot located in a position 
close to a straight line shows that the model sufficiently 
for TOC removal by UV/H2O2 process. The experimen- 
tal and predicted values for Y are shown in Figure 2. Ob- 
servations indicate a very good correlation between the 
results obtained by experiments and the values pre- 
dicted by the statistical model, which shows the success 
of this model. 

Effect of independent variables on TOC removal 

According to Table 3, the Unear effects of four inde- 
pendent variables (initial concentration of TOC, initial 
hydrogen peroxide concentration, pH of solution and 
processing time) are significant. Thus, each variable in 
turn can affect the TOC removal efficiency. P value also 
shows that the effect of pH and initial concentration of 
TOC (P <0.0001) is greater than time and H2O2 concen- 
tration (P <0.001). Interaction and second-degree effects 
of pH and initial TOC concentration are significant, 
confirming the influence of these two parameters. The 
effect of the independent variables on TOC removal 



efficiency as a three-dimensional response surface plot 
is shown in Figure 3. From Figure 3a and c, by decreasing 
the initial concentration of TOC (from 10 to 4 mg/L) the 
removal efficiency increases. Maximum TOC removal was 
achieved at the lowest initial TOC concentration of 4 mg/L. 
Generally, UV radiation and hydrogen peroxide produces 
free radicals such as hydroxyl radicals [1,26,27]. Increasing 
in initial concentration of humic substances generally 
decreases the process efficiency, probably because of 
competition between metabolites of humic substances 
to react with the hydroxyl radicals as a non-selective 
agent. Therefore, increasing the initial concentration of 
humic substances, resulting in increased concentrations 
of intermediate products in the environment, conse- 
quently a significant amount of hydroxyl radicals for the 
degradation of intermediate products is used [27,28]. 
NOM can also compete with H2O2 for reacting with UV 
radiation. Thus, low concentrations of NOMs allow more 
H2O2 photolysis and consequently more hydroxyl radicals 
production [26]. As it is shown in Figure 3a and d, TOC 
removal efficiency increased with increasing H2O2 con- 
centration, in the concentrations range of H2O2 (100 to 
180 mg/L), at low concentration of TOC (4 mg/L), with 
increasing concentrations of hydrogen peroxide, removal 
efficiency increased from about 74% to 82%. Oxidation 
rate increases with increasing H2O2 concentrations. 
This is due to the fact that high concentrations of H2O2 
and more absorption of UV rays, increases hydroxyl 
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Design-Expert® Software 
TOC Removal 

Color points by value of 
TOC Removal: 



130 



95 

90 4 

80 
70 



Normal Plot of Residuals 





J"- 



-0.86 -0.11 0.63 

Internally Studentlzed Residuals 



Figure 1 Normal probability plot TOC removal efficiency (Y). 



radical production [27]. Based on the results of similar 
studies, with increasing concentrations of hydrogen 
peroxide, oxidation efficiency increases to a certain ex- 
tent, but in case of excessive concentrations of H2O2, 
little effect and even in some cases loss oxidation rate 
have been reported. Probably high concentrations of 
hydrogen peroxide act as radical scavenging agent and 
reduce the amount of active radicals [13,29]. Further- 
more, an increase in H2O2 concentration results in its 
direct reaction with hydroxyl radicals and HO2 produc- 
tion with lower oxidation capacity than OH° [27]. Al- 
though the concentration of hydrogen peroxide plays an 
important role in the oxidation process, but in high con- 
centrations due to scavenging the hydroxyl radicals, it is 



less effective [5]. TOC removal efficiency increases with 
the increase of solution pH from 3 to 7 (Figure 3b, c, d). 
And then from pH, 7 to 11 TOC removal efficiency 
decreases. Therefore, the maximum TOC removal effi- 
ciency was achieved at neutral pH. Initial pH of a solu- 
tion plays an important role in the oxidation process. 
Hydrogen peroxide is stable in the pH range of 5-9, the 
more free radicals (H02-, OH-) can be formed by UV 
photolysis. Presence of bicarbonate and carbonate alka- 
linity (especially carbonate alkalinity) can compete with 
NOM for scavenging produced hydroxyl radicals, result- 
ing in decrease of NOM removal at alkaline pH values 
[1,10,26]. The effect of pH on UV/H2O2 process de- 
pends on the nature of the contaminants. For example 



Design-Expert® Software 
TOC Removal 

Color points by value of 
TOC Removal: 

184 
30 



Predicted vs. Actual 




Figure 2 Comparison between the experimental values and the predicted values of RSM model, TOC removal efficiency (Y). 
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Figure 3 Three-dimensional response surface plots for TOC removal versus independent variables, (a) H2O2 concentration and initial TOC 
(b) pH and time (c) Initial TOC and pH (d) H2O2 concentration and pH. 



it has been reported that the decomposition of phenol 
by photochemical process is more effective at acidic pH 
values [30,31]. 

With increasing processing time from 10 to 30 min, 
removal of TOC increased and the maximum removal 
efficiency was obtained at 30 min (Figure 3b). In fact, with 
increasing reaction time, hydrogen peroxide has more 
time to react with UV, resulting in the production of more 
hydroxyl radicals [7,27,28]. It was revealed that the effects 
of pH and initial TOC concentration were more signifi- 
cant than H2O2 concentration and processing time. 

Process optimization and evaluation 

For process optimization, the maximum TOC removal 
efficiency was optimized at minimum consumption of 
H2O2. Optimum values of process variables and test results 
under optimal condition are presented in Table 4. Accord- 
ingly, TOC removal efficiency of 78 .02% was determined 
at the optimum conditions of independent variables. The 
maximum TOC removal efficiency using confirmation 



experiments was 76.50%. Therefore, there is good agree- 
ment between the predicted and experimental affirmation 
results at optimum condition that confirms the developed 
model. It is concluded that RSM is a powerful tool to de- 
termine the exact values of the independent variables [32]. 

Treatment raw water samples under optimal conditions 

In this phase of the study, the TOC removal efficiency 
using UV/H2O2 process for natural raw water samples 
was conducted under optimal conditions. Based on the 

Table 4 Optimum value of the process variables for 
maximum TOC removal efficiency and Predicted and 
experimental value for the responses at optimum 
conditions 



H202 

concentration 
(mg/L) 


pH 


Time 
(min) 


Initial TOC 

(mg/L) 


TOC removal (%) 










Predicted Experimental 


100 


6.12 


22.42 


4 


78.02 76.50 
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values obtained for the optimal concentration of hydro- 
gen peroxide (100 mg/L), time (22.4 m) and normal raw 
water quality (TOC = 4.20 mg/L, pH = 8.15) . The results 
revealed that TOC removal under optimum conditions 
for natural water samples decreased to 62.15% comparison 
with synthetic samples. This difference is probably due 
to the presence of other components (e.g. alkalinity) in 
natural water [1,10,26]. 

Conclusions 

In this work, RSM was employed for optimization of 
TOC removal using UV/H2O2 process. By using the 
Box-Behnken method, four main parameters including 
initial TOC concentration. Hydrogen peroxide concen- 
tration, pH and time were examined. A second-order 
polynomial model was developed using multiple linear 
regression analysis. Statistical test (ANOVA) indicated a 
good agreement between experimental data and the 
built model (R^ = 0.98). The optimal operating conditions 
were determined using numerical optimization tech- 
niques. For this purpose, the maximum TOC removal was 
optimized for minimum consumed H2O2. Accordingly, 
the removal efficiency of TOC, 78.02% in terms of the in- 
dependent variables was optimized. Further confirmation 
experiments under optimum operating conditions showed 
a 76.50% of TOC removal and confirmed that the model 
is accordance with the experimental data. The efficiency 
of TOC removal from natural water based on RSM 
optimum condition was 62.15%. This study showed that 
RSM based on Box-Behnken method is a useful tool for 
optimizing the operating parameters for TOC removal 
using UV/H2O2 process. 
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